Universality of derealization in unconventional dirty superconducting wires with broken 

spin-rotation symmetry 
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In dirty superconducting wires, quasiparticle states at the Fermi level need not be exponentially localized if 
spin-rotation symmetry is broken [Brouwer et al, Phys. Rev. Lett. 85, 1064 (2000)]. Here we present evidence 
that not-localized states are generic in the thick-wire limit, while for wires of finite thickness derealization 
requires fine tuning of the disorder, consistent with earlier results of Motrunich et al. [Phys. Rev. B 63, 224204 
(2001)]. The thick- wire limit is defined as the simultaneous limit where the length L of the wire and the number 
N of propagating channels at the Fermi energy are both taken to infinity with their ratio held fixed. 

PACS numbers: 74.25.Fy, 72.15.Rn, 73.20.Fz, 73.23.-b 
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Localization properties of weakly disordered normal-metal 
wires depend first and foremost on the fundamental symme- 
tries of the wire: presence or absence of time-reversal symme- 
try and spin-rotation symmetry. In addition to these, two more 
symmetries can play an important role in the characterization 
of quantum wires: chiral symmetry and particle-hole symme- 
try. The former is relevant for lattice models with randomness 
in the hopping amplitudes only, while the latter applies to su- 
perconducting wires or normal-metal wires in the proximity 
of a superconductor. 

Is symmetry the only player determining the localization 
properties? For normal-metal wires, it is generally agreed 
on that this is the case when the disorder is weak. For su- 
perconducting wires, different views have been published in 
the literature. Advocates of field-theoretic, random-matrix, or 
Fokker-Planck approaches have mostly considered the sym- 
metry classification as sufficient, thus arriving at a division of 
superconducting wires into four universality classesi^iii^i&i 7 . 
These are labeled C, CI, D, and Dili, and are characterized by 
the presence or absence of time-reversal symmetry and spin- 
rotation symmetry. 1 For the classes D and Dili, where spin- 
rotation symmetry is broken, it was found that wave functions 
are "critical" if disorder is weak, rather than exponentially 
localized, on length scales beyond a characteristic crossover 
scale Nl, where TV is the number of propagating channels in 
the wire and I is the mean free path. Here "critical" indicates 
algebraic decay of the average conductance, opposed to ex- 
ponential decay if wave functions are localized. Disorder is 
"weak" if kpl 3> 1, kjt being the Fermi wave vector. On the 
other hand, using a numerical transfer matrix method to cal- 
culate localization lengths in single-channel wires and strong- 
disorder renormalization-group arguments, Motrunich et al. 
found that, if spin-rotation symmetry is broken, the localiza- 
tion length only diverges if the disorder is fine tuned. 8 They 
concluded that symmetry alone is not sufficient to character- 
ize disordered wires with particle-hole symmetry. Such a find- 



ing is consistent with the expectation that, in a single-channel 
wire, weak spin-orbit scattering or weak superconducting cor- 
relations cannot destroy the localized phase that exists in the 
presence of spin-rotation symmetry and without superconduc- 
tivity. 



In this paper we take a more detailed look at this appar- 
ent contradiction for the case of class D where, in addition 
to spin-rotation symmetry, time-reversal symmetry is broken 
as well. We consider both regimes L> Nl and L < Nl. For 
thick wires (N ^> 1), the regime / <C L <^Nl corresponds to 
the diffusive regime. For thin wires (N of order unity), one 
has Nl ~ I, and there is no diffusive regime. Our main con- 
clusion — to be supported by numerical and analytical ar- 
guments below — is that, as argued in Ref. 8, for quantum 
wires of finite thickness and broken spin-rotation symmetry, 
fine tuning of the disorder is required in order to obtain a di- 
verging localization length. However, we shall also argue that 
in the limit N — > °° of "thick" quantum wires symmetry alone 
is sufficient to determine the localization properties. 9 We note 
that this criterion is not different from that in normal wires. 10 
Thus, as far as the classification based on symmetry is con- 
cerned, superconducting disordered wires are no less univer- 
sal than ordinary normal-metal wires. In the final analysis, 
realization of the quasi-one-dimensional critical behavior of 
class D demands fine tuning of microscopic parameters for a 
wire of finite thickness but becomes generic in the thick-wire 
limit provided the disorder is sufficiently weak. We reach this 
conclusion by studying the competition between localization 
in the standard unitary class (no superconductivity) and crit- 
icality in class D that occurs upon increasing the strength of 
superconducting correlations in the wire. However, the same 
scenario applies to the crossover between the localized phase 
of class C and the critical phase of class D upon breaking of 
the spin-rotation symmetry. 
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As in Ref. 7, our starting point is the Hamiltonian 

?{ = !K, + V, =G0®y0®t3® *NiVFd x , (1) 



In the special case Ia = l v , i.e., in the pure class D, the en- 
semble average dimensionless conductance (g) in the diffu- 
sive regime L<SiNl readsi 



where Oo is the 2x2 unit matrix in the spin grading, yo is the 
2x2 unit matrix in particle-hole grading, and T3 is the Pauli 
matrix in left/right mover grading. The kinetic energy 3C de- 
scribes the propagation of right and left moving quasiparticles 
in N channels at the Fermi level. The "potential" V (x) is an 
%N x %N matrix that accounts for the presence of both disorder 
and superconducting correlations. In particle/hole (y) grading 
it reads 



V = 



v A 

-A* -v T 



(2) 



where v (A) is a Hermitian (antisymmetric) AN x AN ma- 
trix, representing the impurity potential (superconducting or- 
der parameter). The form of the potential V ensures 
that the Hamiltonian !H obeys particle-hole symmetry: H = 
— yi.W yi»i In addition, 9{ (and hence 1/) may obey time- 
reversal invariance 9{ = 1 ?(*T~ l , with 1 = i%\ ®<52- 

In Ref. |7| it was assumed that the potentials A and v in Eq. 
are Gaussian white-noise potentials with the same vari- 
ance. It was for this choice of disorder that the "critical" con- 
ductance statistics was found. Motrunich et al. have argued 
that this choice is special: As soon as the variances of A and 
V are no longer equal, the conductance will decay exponen- 
tially with length L. In addition, they showed that if a term 
proportional to vf^i/21 s is added to the Hamiltonian, which 
corresponds to a staggering of the hopping amplitude in a lat- 
tice version of the problem, 11 the localization length is finite 
even if v and A have the same variance. 

Below we shall analyze the effects of a different variance of 
A and v and of staggering on the localization length, and dis- 
cuss how these effects change with the thickness of the wire. 
We will first address the diffusive regime, where an analytical 
treatment is possible, and then the localized regime, where we 
support our arguments with numerical simulations. 

Different variances of A and v. For simplicity, we assume 
that time-reversal symmetry and spin-rotational symmetry are 
completely broken. Hence we take the potentials v and A as 
Gaussian white-noise, with vanishing means and with vari- 
ances 
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(A^AiM = ^_(8.,8 y7 -8^)S(*->0. 



Here l v and I a are the mean free paths for scattering from po- 
tential disorder and from (fluctuations of) the superconducting 
pair potential A, respectively. When I a — l v , the exact solution 
of Ref. 7 applies. We will refer to this special case as the pure 
class D and consider the generic case Ia ^ l v as a point in the 
crossover between the standard unitary ensemble and class D. 
(The standard unitary ensemble corresponds to the absence of 
superconducting correlations, Ia — <*>.) The combined mean 
free path I for scattering from either v or A is = Z" 1 +17 . 
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For L ^> Nl, the average conductance decays algebraically, 
whereas the "typical conductance" exp((lng)) decreases pro- 
portional to exp[-(2L/7tM) 1 / 2 ]|i 



2L 
kNI 



= -\ ZT77- (5) 



For the standard unitary ensemble Ia = °° one has (g) 
A XI I ■ o(L/Nl) in the diffusive regime, while in the local- 
ized regime (g) °c e - L / ml , {lag) = -L/2NI + 0(1)* 

Let us now consider the generic case ^ Ia < °°, i.e., a 
point in the crossover between the standard unitary ensemble 
and the pure class D. At a fixed but large Ia ^> l v , this same 
crossover can also take place as a function of the wire length 
L. In that case, short L corresponds to the standard unitary en- 
semble, while class D statistics might be found for large L. To 
find the relevant length scale for that crossover we calculate 
the average conductance in the diffusive regime. Following 
the method of Ref.ll3| we derive scaling equations for the av- 
erages of certain traces of transmission and reflection matrices 
t and r of the wire. For Ia,L 3> I and up to corrections of rel- 
ative order (L/Nl) the relevant scaling equations read 
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where we abbreviated ^a = \/IIa/8 and Y3 is the Pauli matrix 
in particle-hole grading. Solving Eqs. (|6) we find the average 
conductance {g) = (txt^t) for Ia,L ;§> 



(g) = 



ANl 



L 



L 2 



Nl 



(7) 



The length scale ^a marks the crossover between the stan- 
dard unitary class and class D: For L <C £,a we recover the 
result (g) =ANl/L+o(L/Nl) of the standard unitary ensem- 
ble, while for L>^ Eq. (0 simplifies to the class D result 
(|4). As Eq. |7} is only valid in the diffusive regime L <C Nl, 
one can only expect to observe the class D result © if ^a is 
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smaller than Nl. In the thick- wire limit this poses no real con- 
straint, and we conclude that, in the thick-wire limit, any finite 
mean free path / A will eventually give rise to a weak localiza- 
tion correction to the conductance according to class D if the 
wire is sufficiently long. 

To find what happens in the localized regime we now com- 
pare the crossover length scale c^ A derived above and the local- 
ization length 4NI in the standard unitary class. If c^ A ^ Nl, 
the wave functions are still governed by the unitary class on 
the scale Nl, not by class D, when localization sets in. Hence 
there is no room for class D physics to play a role and the lo- 
calization length is that of the standard unitary class. On the 
other hand, if cj A "C Nl, the statistics of wave functions and 
conductance is given by that of class D at scales cj A and be- 
yond. Therefore if cj A -C Nl we expect that the wave function 
and conductance statistics for L 3> Nl is critical, as in class 
D. If Nl and c^ A are comparable, the conductance statistics in 
the localized regime is best characterized as in the "crossover" 
between the two universality classes: The localization length 
is expected to be larger than 4NI, but still finite. 

In order to make these statements more precise, we look at 
the normalized inverse localization length, 



X(o) 



o = 



2NI 



(8) 



In the standard unitary ensemble, one has % — > 1 as L — » °° 
with N, Ia, and l v held fixed; critical conductance statistics 
implies %(°°) = 0. A finite value < %(°°) < 1 corresponds 
to exponentially localized wave functions, but with a localiza- 
tion length that is larger than in the standard unitary ensemble. 
Taking the limit L — » °° at fixed values of N, Ia, and l v cor- 
responds to a point somewhere in the crossover between the 
standard unitary ensemble and class D if I a ^ l v . As such a 
choice of parameters corresponds to a. finite distance from the 
pure class D, we expect x(°°) > in this case, corresponding 
to exponentially localized wave functions. This agrees with 
the claim of Motrunich et a/., 8 who found critical statistics at 
the point I a = l v only. 

In the thick-wire limit, however, the picture is different. 
The thick-wire limit is defined as the simultaneous limit L — > 
°° and N — > °°, keeping o = L/2NI fixed. The limit N — > °° 
arises as a semiclassical limit, when the Fermi wavelength is 
much smaller than the diameter of the wire, or by increas- 
ing wire diameter. In the latter case, the diameter should 
not exceed the transverse localization length. Based on our 
calculation for the diffusive regime, we surmise that in the 
thick-wire limit the crossover between the standard unitary 
class and class D is described by a two-parameter scaling 
function %*(c,(X), which is obtained from % by taking the 
limit L,N — > °°, at fixed values of the scaling variables rj and 
a = 2jV//^a- 15 The second scaling variable a is the ratio of 
the localization length 4NI in the standard unitary class and 
twice the diffusive crossover length scale cj A (the proportion- 
ality constant is arbitrary as long as it is of order 1). As a func- 
tion of a, %* interpolates between the standard unitary ensem- 
ble for a = and the critical conductance statistics of class D 
if a — > oo. Note that taking the thick- wire limit at any finite 
and fixed ratio of c^ A and I (or, alternatively, of I a and I) cor- 
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FIG. 1 : Normalized inverse localization length % = — (lng)/o versus 
a = L/2NI. Solid curves are for l A = 201 / (£ A = 5.01 1) and N = 4,6, 
from top to bottom. Dotted and dashed curves are for / A = 101/ 
(£ A = 3.6/) with N = 2,3,4,6,8, and / A = 26/ (£ A = 1.8/) with N = 
2,3,4,6, respectively. Thick dash-dot curves at the top and bottom 
are for the standard unitary class (Z^ 1 = 0) and N = 1 and for the 
pure class D (/ A = 21) and N = 3, respectively. Curves with equal 
values of the scaling variable a = 2M/i; A but different | A /Z or N 
almost coincide; the remaining small differences between them are 
attributed to finite-iV effects. Corresponding values of a are listed to 
the right of each set of curves with the same a. 



responds to the case a = °°, i.e., automatically realizes class 
D. Thus, in the sense of the scaling limit defined above, the 
standard unitary fixed point is "unstable," whereas the class D 
fixed point is "attractive." It is only by allowing cj A // to scale 
at least as fast as N that one can probe non-class-D physics in 
the thick-wire limit. 

We verified the above scenario with numerical simulations 
of % as a function of o for a discretized version of the model 
0. The discretized model is obtained from Eq. by the 
replacement of the continuous potential 1? (x) by the sum 
£ j Vj8(x — a j), where a is a microscopic length scale. The 
delta functions 8(x — y) in Eq. l|3} are replaced by Kronecker 
deltas a 8jf. In Fig.^ results are shown for different values 
of I a/ I and for different N. We took an average over > 2000 
realizations of the disorder potential. Our conclusions are con- 
firmed by the numerical data: First, curves for the same value 
of a = 2M/^ A = 4N(2/// A ) 1 / 2 approach a well-defined limit 
for large N, confirming that a = 2NI/£,a is indeed the true 
crossover parameter. Second, in the crossover between the 
two classes, one still observes localized behavior, but with a 
significantly enhanced localization length, which increases as 
a is increased. 

Staggering. Staggering is described by the addition of a 
term vpTi/lls to the Hamiltonian Q. Without the potential 
term V in Eq. 0, the effect of staggering would be to create 
a gap of size Vf / l s around the Fermi energy e = 0. The effect 
of a random potential V is to create states inside this gap. 
Below we address the localization properties of such states at 
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the Fermi level e = if the random potential V is of the form 
with nonzero superconducting correlations A and without 
spin-rotational invariance and time-reveral symmetry. Based 
on a combination of analytical and numerical arguments, we 
argue that eigenstates at the Fermi level remain critical (as 
in the absence of staggering), provided that the staggering is 
sufficiently weak [see Eq. dlOl below] and the thick-wire limit 



Although the scaling equations l|9} do not form a closed set, a 
solution can be constructed as a series expansion in l/l s . For 
sufficiently weak staggering, 

h > I, (10) 

and for lengths L>/, the perturbation series inl/l s truncates 
after the second order, and one finds the simple solution 

(g) = <tr(l-K M )) = 1 + 0(1), (11) 
where a is now defined as 

L 

a= 2M[i-(///,.) 2 r (12) 

Hence, in the thick-wire limit, staggering causes a renormal- 
ization of the scaling length a, but does not remove the exis- 
tence of a diffusive regime. If staggering exceeds the critical 
strength (1101 . there is no diffusive regime in which (g) is in- 
versely proportional to wire length at e = 0. The existence of 
a critical staggering strength is consistent with numerical sim- 
ulations of periodic-on-average single-channel quantum wires 
by Deych et al~ 

Based on this result for the diffusive regime, we propose 
that in the simultaneous limit N — » °°, L — > °° the normal- 
ized inverse localization length % = — (lng)/o approaches the 
same scaling limit %* (a) as in the absence of staggering, if a 
is given by Eq. (I12> . In Fig.|2]we show the result of numerical 
simulations for the discretized version of Eq. Q with stagger- 
ing. We have compared the cases I v /Ia = 1 or 0, with or with- 
out (weak) staggering and find that for sufficiently large N all 
curves % versus o with the same symmetry of the Hamiltonian 
collapse. That is, curves with I v /Ia = 1 (or, in fact, any fi- 
nite ratio I v /Ia) approach the class D result, while curves with 
Iv/Ia — approach the result of the standard unitary class. 
Small N curves are slightly different, which again indicates 
the absence of universality when a diffusive regime is absent. 



is taken. 

To study the effect of staggering on the e = eigenstates 
of the Hamiltonian Q, we first investigate the conductance 
in the diffusive regime. Following the method of Ref. 13|, we 
construct scaling equations for the moments R nM1 = (tr r n r^ m ). 
To leading order in N, they read 



(9) 



In conclusion, our analytical and numerical results confirm 
the claim of Motrunich et al. that the conductance statistics 
in the localized regime is not determined by symmetry alone 
for a superconducting quantum wire with a small number N 
of propagating channels at the Fermi level. 8 Such universal- 
ity is only achieved for "thick" quantum wires with N 3> 1. 
Which situation is relevant depends on the physical system at 
hand. For quasiparticle states at the edges of superconducting 
planes in high-r f superconductors, the nonuniversal scenario 
of small N may be more applicable, while we expect that the 
universal scenario of large N is more appropriate for dirty un- 
conventional superconducting wires or normal-metal wires in 
the proximity of a superconductor. Possible relevant pertur- 
bations to the critical behavior in the thick-wire limit of the 
superconducting class D are disorder sufficiently strong to in- 
duce transverse localization, restoration of the spin-rotation 
symmetry, or breaking of the particle-hole symmetry (e.g., by 
tuning the chemical potential away from the the Fermi level 
or by electron-electron interactions). All of these perturba- 
tions are expected to restore the localized phase at large length 
scales. 

We would like to point out that the universality — or lack 
thereof — for the localization properties in superconducting 
quantum wires is not different from the case of normal-metal 
wires. In both cases, universality is expected in the thick- 
wire limit only. Or, in the language of scaling flow: In both 
cases, the fixed point with less symmetry is attractive, but the 
flow towards the fixed point is only complete if the thick-wire 
limit is taken. In that respect, note that for the crossover be- 
tween the standard unitary class and the superconducting class 
D considered in this paper, the superconducting class is the 
class of less symmetry. That fact constitutes a qualitative dif- 
ference between class D and the three chiral classes, which 
also may exhibit critical conductance statistics for wire length 
L 3> M»ii In contrast to the superconducting class D, the chi- 
ral classes have very high symmetry, which makes them rather 
unstable and which will render observation of the critical con- 
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ductance statistics in these classes more difficult. 



FIG. 2 : Normalized inverse localization length % = — (In g) /a versus 
scaled wire length a with and without staggering. With staggering, 
o is given by Eq. <12> . without staggering o = L/2NI. Upper panel 
shows results for normal quantum wires (l\ = °°). Solid curves are 
for l v = l = 0.80/s and for N — 4, 6; dashed curves are without stag- 
gering and N = 4, 6. The slight difference (less than 5%) between 
curves with and without staggering is attributed to a small deviation 
from Eq. <12> due to the finite disorder strength in the simulations. 
As all curves for different N collapse, no single curves are labeled. 
Lower panel shows results for superconducting quantum wires with 
/a = lv Solid curves are for l v = l\ = 21 = 1.60/ s and N = 1,2,3,4, 
and 6. Dashed curves are without staggering and N = 3, N = 4. They 
almost coincide with the curve with staggering and N = 6 and are not 
labeled individually. 
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